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' Abstract. Phase- locking governs the phase noise in classical clocks through effects de- 

04 , scribed in precise mathematical terms. We seek here a quantum counterpart of these effects 

(~| ' by working in a hnite Hilbert space. We use a coprimality condition to define phase-locked 

^ , quantum states and the corresponding Pegg-Barnett type phase operator. Cyclotomic sym- 

metries in matrix elements are revealed and related to Ramanujan sums in the theory of 
prime numbers. The employed mathematical procedures also emphasize the isomorphism 
between algebraic number theory and the theory of quantum entanglement. 

(N 

Time and phase are amongst the most federal concepts of science. Hersh mentions Kant as 
putting time before number and making the intuition of time the origin of arithmetics [1]. 
, The present work can be considered as an extension of a longstanding effort to model phase 

[ noise and phase- locking effects that are found in highly stable classical oscillators. It was un 



' ambiguously demonstrated that the observed variability (i.e., the 1/f frequency noise of such 

oscillators) is related to the finite dynamics of states during the measurement process and to the 
precise filtering rules that involve continued fraction expansions, prime number decomposition, 
and hyperbolic geometry [2]- [4]. We are interested here in possible quantum counterparts of 
^ ■ these effects. The problem of defining quantum phase operators was initiated by Dirac in 1927 

I [5]. After more than three quarters of a century, it is still an interesting and essentially open 

problem. For an excellent review, see [6]. Our starting point is the Pegg and Barnett (PB) 
quantum phase formalism [7]- [8] where the calculations are performed in a Hilbert space Hg of 
finite dimension q [9]. In the PB framework the phase states are the so-called discrete quantum 
^ I Fourier transforms (QFT) of the number states, i.e., the following superpositions 



o 



q-l 



-V.j;expf^)|„>. (1) 



n=0 

The states \9p) form an orthonormal set and in addition the projector over the subspace of phase 
states is J2l=o I^p)(^pI = 1<? where Ig is the identity operator in Hg. The inverse QFT implies 
\n) = q^^l"^ Z]p=o 6xp(— ^^y^)|^p). As the set of number states |n), the set of phase states |0p) is 
a complete set spanning Hg. In addition, the QFT operator is a g x (7 unitary matrix of entries 
4?) = ^exp(2ml|^). 

Prom now on, we restrict ourselves to the class of phase states \Q'^ for which p and q 
satisfy the coprimality condition (p, = 1, where (p, g) is the greatest common divisor of p 
and q. Differently from the phase states (1), form an orthonormal base of a Hilbert space 
whose dimension is given by the number of irreducible fractions "pjq. This number is known in 
Mathematics as the Euler totient function <^{(i). 

^e-mail: planat@lpmo.edu; Fax: 33(0)381853998 
^e-mail: licr@ipicyt.edu.mx; Fax: 524448335412 



1 



Guided by the analogy with the classical situation [2] we call the irreducible states the 
phase-locked quantum states. In the space of quantum states, they generate a cyclotomic lattice 
L [10]-[14] with a generator matrix M of matrix elements Kpn and (p, q) = 1 and size (p{q) x (p{q). 
We notice that the corresponding Gram matrix of L, namely H = M^M, has matrix elements 
^nl = Cg(n — Z), which are Ramanujan sums 

(2 (2, 

(p!«)=i 

where r = q/{q,n). Ramanujan sums are thus defined as sums over the primitive characters 
exp(2i7r^), {p,q) = 1, of the group Zq = Z/qZ. In the equation above, /x(r) is the Mobius 
function which is if the prime number decomposition of r contains a square, 1 if r = 1, 
and (—1)^ if r is the product of k distinct primes [15]. Ramanujan sums are relative integers 
that are quasiperiodic in n with quasi period (p(q) and aperiodic in q with a type of variability 
imposed by the Mobius function (introduced by Ramanujan in the context of the Goldbach 
conjecture). Recently, we demonstrated that Ramanujan's sums are also useful in the context of 
signal processing as an arithmetical alternative to the discrete classical Fourier transform [16]^. 

The projection operator over the subset of phase-locked quantum states \9p) is calculated as 
the coprimality-constrained sum of one-state projectors, i.e., 

piock^ g g |^;)(0;i = lgc,(n-O|n)(/|, (3) 

P=o p=o ^ n,l=0 

(P,9)=l (P,9)=l 

so that the matrix element of the projection is such that g(n|Pg|/) = Cq{n — I). This is similar 
to the elements of the cyclotomic Gram matrix and sheds light on the equivalence between 
cyclotomic lattices of algebraic number theory and the quantum theory of phase-locked states. 
The quantum phase-locking operator is defined 

q-l 

©5°''^= E WK\ = ^Pi^ where = 27r^ . (4) 

p=() ^ 
(p,y)=i 

The Pegg and Barnett phase operator Qq [7]- [8] results from the removal of the coprimality 
condition {p,q) = 1. Using the number operator Nq = X)^=ora|n)(n| a generalization of Dirac's 
commutator [Qq,Nq\ = —i is obtained. 

In contrast, the phase number commutator for phase-locked states as calculated from (4) is 

C^'^ = [0r,^^.] = ^ EO-riK(n-Z)|n)(Z|, (5) 

^ n,l=0 



^In the discrete Fourier transform the signal processing is performed by using all roots of unity of the form 
exp{2inp/q) with p from 1 to g and taking their n**" power ep(n) as basis function. We generalized the classical 
Fourier analysis by using Ramanujan sums Cq{n) as in (2) instead of ep(n). This type of signal processing is more 
appropriate for arithmetical functions than is the ordinary discrete Fourier transform, while still preserving the 
metric and orthogonal properties of the latter. Notable results relating arithmetical functions to each other can 
be obtained using Ramanujan sums expansions while the discrete Fourier transform would lead instead to low 
frequency tails in the power spectrum. 
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and has antisymmetric matrix elements (Z|Cg°^'^|ra) = ^{l — n)cq{n — I). 

The finite and cyclotomic quantum mechanical rules are mostly encoded in the expectation 
values of the phase operator and phase variance (see Figures 1-3). 
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Fig. 1 

The phase expectation value versus the dimension q of the Hilbert space. Plain lines: /? = 0. Dotted 

lines: (3=1 
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Fig. 2 

The phase expectation value versus the phase parameter /J. Plain lines: q = 15. Dotted lines: q = 13. 
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Fig. 3 

The phase variance versus the dimension q of the Hilbert space. Plain hnes: /3 = 1. Dotted hnes: 
(3 = TV. The horizontal line at w^/S is the variance of the classical phase. 

Rephrasing Pegg and Barnett, we consider a pure phase state |/) = J2^^oUn\n) having u„ 
of the form ^ 

Un = — exp(in/3), (6) 

where /3 is a real phase parameter. Defining the phase probabiUty distribution as {Op\f)'^; the 
phase expectation value as (0)^°'^'^) = J2p(^p{^p\f)'^'i ^^'^ the phase variance as {AQ^y^^^ = 
^p{Op — {Q^q'^^))'^ {0'p\f)'^ , where the summations are taken from p = to g — 1 and {p,q) = 1; 
the following result is obtained 

i^T'') = ^ - n) exp[i^(n - I)], (7) 

(^@2)lock ^ 4^@lock^ + l^((eiock^ _ 27r), (8) 

with the modified expectation value 

9-1 

(eiock^ = _ ^ - n) exp[i/3(n - /)], (9) 

^ n,l=0 

and the modified Ramanujan sums 

Cq{m) = ^ {p/qf exp (2iTrm-) . (10) 

(p,«)=i 

Fig. 1 illustrates the phase expectation value versus the dimension q for two different values 
of the phase parameter /?. For /? = 1, there are peaks at dimensions q = p^, which are powers 
of a prime number p. For /? = 0, the peaks are smoothed out due to the averaging over the 



4 



Ramanujan sums. Fig. 2 shows the phase expectation value versus the phase parameter (3 for 
two dimensions of the Hilbcrt space, tj' = 15 and q = 13, respectively. For the case of the prime 
number q = 13, the phase expectation value is of considerable higher level than for q = 15, with 
absorption-like lines at isolated values of (3. For g = 15, which is not a prime power, the phase 
expectation value is at a much lower level and more random than for the case of the neighbour 
prime number. 

Fig. 3 illustrates the phase variance versus the dimension q. The case of /3 = 1 also leads to 
peaks at prime powers. Similar to the expectation value in Fig. 1, this case is thus reminiscent of 
the Mangoldt function A(n) defined as Inp if n is the power of a prime number p and otherwise. 
This function arises in the frame of prime number theory [2] in the formula for the logarithmic 
derivative of the Ricmann zeta function = X^n^o mean value oscillates around 1 

with an error term which is explicitely related to the positions of Riemann zeros of Q{s) on the 
critical line s = \. The error term shows a power spectral density close to that of 1// noise [2] . 
Finally, the phase variance is considerably smoothed out for f3 = it and is much lower than the 
classical limit 7r^/3. Thus, the parameter (3 is useful in defining quantum phase- locked states 
with small phase variances for a whole range of dimensions. 

In conclusion, several interesting properties are pinpointed here as a result of introducing 
cyclotomic discreteness in the Pegg-Barnett quantum phase formalism. The interplay between 
results associated to prime number theory and quantum phase-locking phenomena looks quite 
stimulating. We recall that the idea of quantum teleportation was initially formulated by Bennett 
et al. in finite-dimensional Hilbert space [17]. Yet independently of this, one can conjecture that 
cyclotomic locking could play an important role in many fundamental tests of quantum mechan- 
ics related to quantum entanglement. As a matter of fact, Munro and Milburn [18] have already 
conjectured that the best way to see the quantum nature of correlations in entangled states is 
through the measurement of the observable, canonically conjugate of the photon number, i.e., 
the quantum phase. Moreover, in their paper dealing with the Greenberger-Horne-Zeilinger 
quantum correlations, Munro and Milburn presented a homodync scheme requiring discrete 
phase measurement, yet to be experimentally realized in the ultrahigh detector efficiency limit. 
Nevertheless, a homodyne scheme automatically implies phase locking measurements with cy- 
clotomy playing a role. We also recall that a universal algorithm for the optimal quantum state 
estimation of an arbitrary finite dimensional system in terms of a sequence of projectors has 
been introduced by Derka et al [19]. In one of their two examples, the sequence of projectors is 
very similar to the Pegg-Barnett projectors (in fact they noticed the analogy with the PB phase 
operator). Therefore, the quality of the cyclotomic phase estimation through a direct gener- 
alization of the maximal mean fidelity introduced in their paper could be utilized as well. It 
may be possible that the type of quantum phase-locking described in this paper may have some 
applications to the synchronization of remote oscillators by quantum entanglement as already 
proposed in [20]. 
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